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ABSTRACT 

One of the most essential but uncertain processes for producing cosmic-rays (CRs) and their spectra 
is how accelerated particles escape into the interstellar space. We propose that the CR electron 
spectra at > TeV energy can provide a direct probe of the CR escape complementary to the CR 
nuclei and gamma-rays. We calculate the electron spectra from a young pulsar embedded in the 
supernova remnant (SNR), like Vela, taking into account the energy-dependent CR escape. SNRs 
would accelerate and hence confine particles with energy up to 10 15 5 eV. Only energetic particles 
can escape first, while the lower energy particles are confined and released later. Then the observed 
electron spectrum should have a low energy cutoff whose position marks the age of the pulsar/SNR. 
The low energy cutoff is observable in the > TeV energy window, where other contaminating sources 
are expected to be few due to the fast cooling of electrons. The spectrum looks similar to a dark 
matter annihilation line if the low energy cutoff is close to the high energy intrinsic or cooling break. 
The future experiments such as CALET and CTA are capable of directly detecting the CR escape 
features toward revealing the origin of CRs. 

Subject headings: acceleration of particles - cosmic rays - pulsars: general - supernova remnants 



1. INTRODUCTION 

The origin of cosmic-rays (CRs) is a long standing 
problem since its discovery. The number spectrum of 
the nuclear component of CRs can be fitted with a 
broken power-law: N(e) oc £~ s with s c± 2.7 below 
e ~ 10 15 5 eV (the "knee" energy). The charged CR 
particles are considered to propagate diffusively with an 
energy-dependent manner (the higher energy particles 
can diffuse faster), so the intrinsic CR spectrum at the 
source should be harder than that observed at the Earth. 
Although there are some ambiguities in the energy de- 
pendence on the diffusion coefficient, the source spectral 
index is considered to be around s ~ 2.2—2 .4 (e.g. Strong 
& Moskalenko 1998). 

It is often argued that the CRs with energy smaller 
than the knee energy (or even up to 10 eV) are 
originated in our Galaxy. The most widely accepted 
paradigm for the galactic CR production process is the 
diffusive shock acceleration (DSA) at supernova rem- 
nants (SNRs) . The theory of DSA (for reviews see Bland- 
ford & Eichler 1987; Malkov & Drury 2001) can naturally 
derive the power-law spectrum of particles accelerated in 
the SNR shock. In fact, H.E.S.S. detected TeV gamma- 
ray emissions from the shell of young SNRs (Aharonian 
et al. 2004, 2005), and Fermi and AGILE have detected 
GeV gamma-ray emissions from middle-aged SNRs in- 
teracting with nearby molecular clouds, which are likely 
to be generated via hadronic process, i.e. inelastic colli- 
sions between CR protons accelerated at SNRs and am- 
bient protons (Abdo et al. 2009b, 2010a, b, c; Tavani et 
al. 2010; Giuliani et al. 2010). However, in the conven- 
tional DSA theory the predicted spectral index of CRs 
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accelerated at the shock with a large Mach number is 
s = 2 and harder than that expected from the observa- 
tions. Moreover, the non-linear theories of DSA predict 
a harder source spectrum s < 2.0, which seems to make 
the contradiction between theories and observations even 
worse. 

It has been proposed that these observational facts 
can be interpreted by the model of the energy-dependent 
CR escape from the SNR (Ptuskin & Zirakashvili 2005; 
Drury et al. 2009; Reville et al. 2009; Caprioli et al. 
2009; Gabici et al. 2009; Ohira et al. 2010a; Fujita et al. 
2010; Casanova et al. 2010). In the escape model it is as- 
sumed that the most energetic particles leave the SNR at 
the beginning of the Sedov phase and, as the shock slows 
down and the magnetic field decays, the lower energy 
particles, which have been confined around the shock 
front by the magnetic field, can leave the shock grad- 
ually. This process can largely affect the CR spectrum 
detected at the Earth (Ohira et al. 2010a; Caprioli et al. 
2010). In fact, the spectral index of the CR escaping the 
SNR s csc is different from that of the CR inside the SNR 
s, 

s esc = s + — , (1) 
a 

where the escape energy e osc and the normalization factor 
of the CR production rate are assumed to be proportional 
to t~ a and t@, respectively. For example, in the phe- 
nomenological model by Gabici et al. (2009) a ~ 2.6 (see 
also Ohira et al. 2010a). Therefore the observed spectral 
index can get softer than the conventional value (s ~ 2) 
if the normalization factor of the spectrum gets larger 
with time (i.e. {3 > 0; Ohira et al. 2010a), so we may in- 
terpret the observed CR spectrum in the context of DSA 
with the escape model. Moreover, recent gamma-ray ob- 
servations of middle-aged SNRs interacting with molecu- 
lar clouds by Fermi and AGILE show that those spectra 
are fitted by a broken power-law with a break energy of 
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~ 1 — lOGeV (see the references shown above), and they 
can be interpreted by the energy-dependent CR escape 
from an SNR (Aharonian & Atoyan 1996; Gabici et al. 
2009; Ohira et al. 2010bfl. However, these gamma-rays 
are the secondary emissions of CRs from an SNR, and so 
these observations are only the indirect evidences of the 
CR escape scenario. 

On the other hand, the direct observations of CR 
electrons/positrons have been also greatly advanced. 
The PAMELA satellite discovered the excess of the CR 
positron fraction (Adriani et al. 2009) and, as for the 
flux of CR electrons plus positrons the experiments such 
as ATIC/PPB-BETS (Chang et al. 2008; Torii et al. 
2008b), Fermi (Abdo et al. 2009a; Ackermann et al. 
2010) and H.E.S.S. (Aharonian et al. 2008, 2009) have 
revealed an excess from the conventional model (see also 
Meyer et al. 2010 for the energy calibration between 
Fermi and H.E.S.S.). These results suggest that there 
are some additional electron/positron sources. Possible 
candidates include astrophysical sources such as pulsars 
(Shen 1970; Atoyan et al. 1995; Chi et al. 1996; Zhang 
& Cheng 2001; Grimani 2007; Kobayashi et al. 2004; 
Buesching et al. 2008; Hooper et al. 2009; Yuksel et al. 
2009; Profumo 2008; Malyshev et al. 2009; Grasso et al. 
2009; Kawanaka et al. 2010; Heyl et al. 2010; Blasi k 
Amato 2010), supernova remnants (Shaviv et al. 2009; 
Blasi 2009; Blasi & Serpico 2009; Fujita et al. 2009; Hu et 
al. 2009; Biermann et al. 2009; Mertsch & Sarkar 2009; 
Ahlers et al. 2009), gamma-ray bursts (GRB; Ioka 2010; 
Kistler & Yuksel 2009), microquasars (Heinz & Sunyaev 
2002) and white dwarf pulsars (Kashiyama et al. 2010). 
Dark matter annihilations/decays (e.g. Hooper 2009) 
and the propagation effect (Delahaye et al. 2008; Cowsik 
& Burch 2009; Stawarz et al. 2010) are also the possi- 
ble processes for making the CR electron/positron excess 
(for the comprehensive review, see Fan et al. 2010). 

In the near future within a few years, Alpha Mag- 
netic Spectrometer - 02 (AMS-02) experiment will mea- 
sure the positron fraction up to ~ ITcV (Beischer et al. 
2009; Pato et al. 2010; Pochon 2010) and CALorimetric 
Electron Telescope (CALET) experiment will explore the 
electron spectrurr0 up to ~ lOTeV with an energy reso- 
lution better than a few percent (Torii et al. 2008a). In 
addition, the future Cherenkov Telescope Array (CTA) 
will be able to measure the CR electron spectrum up to 
~ 15TeV (CTA consortium 2010). These experiments 
will open the window to the CR astrophysics with the 
TeV electron and positron components. Especially, as 
Kobayashi et al. (2004) have pointed out, in the TeV 
energy band a few nearby sources may leave spectral sig- 
natures and we will be able to see a spectral shape of 
CR electrons/positrons from a single source, while in the 
lower energy band we can see only the superposed spec- 
trum from multiple sources. 

In this paper, we suggest a possibility that the precise 
measurement of the CR electron spectrum in the very 
high energy band (> 1 - lOTeV, expected with CALET 

3 Recently some authors have tried to explain this spectral break 
by considering the other processes such as the propagation of the 
SNR shock in the molecular clouds (Uchiyama et al. 2010), or the 
two-step acceleration of cosmic-ray particles by the shocks gener- 
ated by the turbulence behind the SNR shock (Inoue et al. 2010). 

4 Hereafter we express the spectrum of electrons plus positrons 
as just "electron spectrum". 



and CTA) could directly prove the energy-dependent CR 
escape if those electrons/positrons are generated by a 
pulsar embedded in an SNR. This is quite a natural sit- 
uation to be realized because a pulsar should be gener- 
ated in the center of a core-collapse supernova and CR 
electrons / positrons are considered to be generated in the 
pulsar wind nebula formed inside the SNR. If the CR es- 
cape really occurs in a young nearby pulsar/SNR system, 
the electron spectrum from it will show a unique feature 
as explained in the following and, when it is observed, 
that will be the first direct evidence of the CR escape 



scenario. 



2. A SIMPLE EXAMPLE 



To illustrate our main idea, we show a clear example 
of the effect of the energy-dependent CR escape in the 
direct electron observations. In the calculations of CR 
electrons/positrons, we usually assume that they are in- 
jected into the interstellar matter (ISM) from a pulsar 
with a spectral form of 



Qe{Se) — Qo^e "i 



(2) 



where e e is the energy of electrons/positrons. The ob- 
served electron spectrum f(e e ,r,t) can be obtained by 
solving the diffusion equation 

-f = D(e e )V 2 f + — [P( e *)f\ + Q e (e e ,r,t), (3) 

where D(e e ) = Dq(1 + e e /3GeV) 5 is the diffusion coef- 
ficient and P(s e ) is the energy loss rate. Here we adopt 
D a = 5.8 x 10 28 cm 2 s _1 , 5 = 1/3 that is consistent with 
the boron to carbon ratio according to the latest GAL- 
PROP code, and P(e e ) = -be\ with b = lO^GeV^s" 1 
which includes the energy loss via synchrotron emission 
and inverse Compton scatterings (with Thomson approx- 
imation). Then, if electrons/positrons are injected from 
a point-like source instantaneously (i.e. Q e (e e ,r,t) oc 
S(t)S(r)), the observed spectrum is simply written as 



Qo£ e 



_3/2 J3 
7T U diff 



(1 - bte, 



\a-2 e -(r/d diS y 



(4) 



where ddiff ~ ^/4D(e e )t is the diffusion length of CR 
electrons/positrons. This spectrum is roughly propor- 

exp(— r 2 / (AD(e e )t)) up to the sharp 



tional to e, 



-a-35/2 



cutoff at e e ~ 1/ (bt) , and exponentially damps beyond 
the diffusion length ddiB (Atoyan et al. 1995; Ioka 2010). 

However, if the energy-dependent escape of CR parti- 
cles from the shock is taken into account, the electron 
spectrum would have a sharp cutoff in the low energy 
side because lower energy CRs cannot escape into the 
ISM. The energy of particles which are marginally capa- 
ble of escaping to the ISM e csc is generally determined 
by the confinement condition of CR particles, i.e. the 
equality between the diffusion length of particles and the 
characteristic size of the system: 



j ^sh(^esc) D 
'cliff = ~ -Ksh, 



(5) 



where Z? s h, w s h and R s h are the diffusion coefficient 
around the supernova remnant shock, the shock veloc- 
ity, and the size of the system, respectively. Then the 



Escaping CR Electron Spectrum 



3 



> 

CD 

(3 



CO 100 r 



CM 

CD 

o 



X 
13 



10 




with escape effect 
w/o escape effect 
ATIC 
H.E.S.S(08) 
PPB-BETS 
Fermi 
H.E.S.S.(09) 



10 1 10 2 10 3 10 4 10 5 

Energy e e [GeV] 

Fig. 1. — The electron spectrum predicted from a transient source 
with the escape energy of £ cac = 7TeV (solid line) added with the 
background model (dotted line), compared with the ATIC/PPB- 
BETS/H.E.S.S. /Fermi data. The spectrum without assuming the 
energy-dependent escape is also shown (dashed line). We assume 
that a source at r = 700pc from the Earth produces pairs a 
time tagc = 1.0 X 10 4 years ago with total energy E e + = E e - = 
0.5 X 10 49 erg and spectral index a = 2.4. 

energy e osc is generally a function of time t due to the 
evolution of Z? s h, R s h and u s h- 

Fig. 1 shows the electron spectrum assuming that 
only electrons and positrons above the escape energy 
e csc are injected instantaneously. The energy of elec- 
trons/positrons which is initially e csc becomes e csc /(l + 
bs esc t aee ) as a result of radiative cooling (where i a ge is the 
time since the emission), and we can clearly see a sharp 
cutoff of the spectrum at that energy, which is the ef- 
fect of energy-dependent escape of CR particles. If such 
a spectrum is confirmed by the future experiments, it 
would be the strong support for the CR escape scenario 
of the particle acceleration process at the SNR, and we 
may also get the information of e C sc at the time of the 
source age i age . 

3. MODEL OF CR ELECTRON /POSITRON ESCAPE 

Let us consider the more sophisticated model than dis- 
cussed in the previous section. We consider a pulsar 
emitting electrons and positrons embedded in the SNR 
(i.e. it has not been evacuated from the SNR by the 
natal kick). A pulsar is considered to be an efficient 
factory, because its rotating magnetic field would pro- 
duce a strong electric field around the pulsar and then 
a large number of e pairs would be produced via elec- 
tromagnetic cascades. The created pairs would stream 
away by a centrifugal force as a pulsar wind that ends 
with a termination shock where the acceleration of elec- 
trons/positrons may occur. Hereafter we assume the 
production rate per energy from a pulsar having a spec- 
trum with a cutoff power-law shape: 



iV e , pr (£ e ,t) = Qo(t)e e a exp 



S-e.cut 



(6) 



where the high energy break is fixed as e e!CUt = lOTeV. 
According to the gamma-ray observations of PWNe, the 
electrons/positrons accelerated up to the energy of 10- 
lOOTeV seem to exist in the nebulae (e.g., Aharonian et 
al. 2006), so this assumption is reasonable. Here Qo(t) 
is given to satisfy 



/CO 
d£ e s e N etPI (£ e , t), 



(?) 



where L ejPr (i) is the electron/positron production lumi- 
nosity and £ e!m i n is set to be lGeV. We assume that this 
luminosity is proportional to the spin-down luminosity 
of the pulsar: 



^ / e,pr(^) <X 



1 



(i + t/T y 



(8) 



where To ~ 10 2 4 year is the spin-down timescale, which 
is related to the surface magnetic field of the pulsar 
(Shapiro & Teukolsky 1983). 

In order to evaluate the CR spectrum in the escape 
scenario, we should assume the time evolution of £ C sc- 
In this study we adopt two models for the functional 
form of £ e sc- In the first model, we assume for £ eS c{t) a 
power-law behavior, and determine its normalization and 
power-law index according to the hypothesis that SNRs 
are responsible for the observed CRs with the energy 
from ~ lGeV up to the knee energy (~ 10 15 5 eV). Then 
e esc (t) should reach the knee energy at the end of the 
free expansion phase (i.e. the beginning of the Sedov 
phase; isedov) and should decrease down to lGeV at t ~ 
10 5 / 2 isodov (i.e. the end of the SNR expansion; Gabici 
et al. 2009; Ohira et al. 2010a): 



,(t) = 10 b - 5 GeV x 



t 



^Scdo 



(9) 



As the second model for the evolution of e csc we 
adopt the one discussed by Ptuskin & Zirakashvili (2005), 
which takes into account the modification of a shock 
structure due to the CR pressure, as well as the non- 
linear dissipation of magnetic turbulence. They solve 
the steady-state equation which determines the energy 
density of the magnetohydrodynamic turbulence W: 

uvw{k) = 2 (r cr (fc) - ri(fc) - r nl (fc)) w(k), (10) 

where u is the flow velocity (here it is equal to the shock 
velocity fi s h), & is the wave number of the turbulence, 
and r cr , T\ and T n i are the wave growth rate at the shock 
due to the CR streaming instability, the damping rate of 
waves in the background plasma due to the ion-neutral 
and electron- ion collisions (linear damping), and due to 
the wave- wave interactions (non-linear damping), respec- 
tively. The mathematical expressions for these functions 
are shown in Ptuskin & Zirakashvili (2005), and by solv- 
ing this equation (while the term of linear damping is 
neglected), we obtain the threshold particle energy for 
escape as a function of the shock velocity w s h. Since we 
know the time dependence of the shock radius and the 
shock velocity in the Sedov phase, we can derive e esc (t) as 
a function of time. When the age of the SNR is younger 
than < 10 5 years, the evolution of e esc (t) in this model 
can be approximated as 



3 (t) ~ 10 5 - 5 GeV x 



^Scdo 



(11) 



Fig. 2 shows the evolutions of e csc (t) in two models 
described above. Once we fix the time dependence of 
e csc , we can describe the number luminosity per energy 
of escaping electrons and positrons as 



-Ne,esc,l( e e, t) = N e>pI (e e , t)B(s e - £ c 



(12) 
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Fig. 2. — The threshold energy for escape from the SNR shock 
£ csc in the Sedov phase as a function of time. The thick solid line 
and the thin solid line correspond to the model by Ptuskin & Zi- 
rakashvili (2005) and Gabici et al. (2009), respectively. The cutoff 
which appears in the spectrum due to the energy loss (including 
the KN effect) during the propagation is also shown (dashed line). 

where &(x) is the step function^]. 

There is another flux component which should be taken 
into account. The electrons/positrons which have the en- 
ergy lower than e csc (t) are confined in the SNR and lose 
their energy adiabatically (Ptuskin & Zirakashvili 2005). 
As e csc (t) decreases with time, some of the confined and 
adiabatically cooled electrons/positrons can escape the 
shock surface when their energy becomes greater than 
£ esc (t). First, the CR electron/positron number per en- 
ergy confined in the SNR can be written as 



N. 



n f(e e , t) — 



dt'N, 



' del 



(13) 



xQ{e esc {t')-e' e ), 

where e' e is the energy of the electrons/positrons at the 
time t' , and it is e e at the time t. 

The adiabatic loss is determined by the expansion law 
of the SNR, 



(fee 

~dt 



-RsNR 
-RsNR 



(14) 



where -Rsnr and -Rsnr are the radius of the SNR shell 
and its expansion velocity, respectively. In the following 
we assume that the SNR is in the Sedov phase, in which 
the time dependence of -Rsnr is expected to be propor- 
tional to i 2 / 5 . Then the energy of electrons/positrons 
confined in the SNR would evolve as 



2/5 



(15) 



and therefore, the distribution function of confined elec- 
trons/positrons can be evaluated as 



N, 



conf 



£Scdo 



x 2/5(l-a) 

dt'Q Q (t')s- a [^ 



x exp 



Se{t/t 



M2/5 



£e.cut 

5 Strictly speaking, the spectrum of the escape flux has a finite 
width around the escape energy e esc . However in the usual case 
the step function is a good approximation (Ptuskin & Zirakashvili 
2005; Caprioli et al. 2009). 
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2/5^ 



(16) 



If the decreasing rate of £ CS c is faster than the adiabatic 
cooling rate, — (2/5)e csc /t, then a part of confined elec- 
trons/positrons can escape the SNR shock. By using the 
same logic in deriving Eq. (21) of Ptuskin & Zirakashvili 
(2005) in the case of expanding media, we can evaluate 
the spectrum of such particles as 



:<2 (e e ,t) = -S(e e 



(£e,t). 



de c 



Of 



2£e 

5t 



xN e , con{ (e e ,t). (17) 

Hereafter we consider the spectrum of escaping elec- 
trons/positrons number luminosity per energy as the sum 
of above two components: 

N e ,csc(£e) = ^Ve,csc,l(£e) + ^Ve,o S c,2 (e e ) ■ (18) 

We neglect the radiative energy loss of elec- 
trons / positrons during the confinement for simplicity. As 
we have shown in Eq. ([TT]) , the electron flux of the second 
component -/V e!CSC .2(£ e , t) is determined by the difference 
between the decline rate of the escape energy £ csc and 
the energy loss rate of electrons/positrons. In the cases 
shown in this study, the decline rate of the escape energy 

is - ae e /t ~3x 10~ 9 (£ e /lTeV) (i/10 4 yr) -1 GeV sec" 1 
(a ~ 1 — 2.6, depending on the model) , while the 
adiabatic cooling rate and the radiative cooling rate is 

- 1.2 x 10~ 9 (£ e /lTcV) (t/lOV)" 1 GeV sec" 1 , - be 2 e ~ 

10 -10 (£ e /lTeV) 2 GeV sec -1 , respectively, where for the 
latter we take the cooling rate for the interstellar space 
(see Sec. 2). Therefore, even if we take into account 
the radiative cooling the flux of the electrons/positrons 
which have once been confined in the SNR does not 
change its order from our calculation. Moreover, in the 
interstellar space the diffusion timescale for TeV elec- 
trons can be estimated as 



tdifl ' 



AD{e e ) 

1.4 x 10 4 yr(-^T 1/3 ( 
y V3TeV/ V 



300pc 



(19) 



and therefore the energy loss of TeV electrons during the 
propagation is at most Ae e /e e ~ 1 — (1 + &<diff£e) 1 ~ 
10%. Then we can say that in this energy range both 
of the energy losses in the PWN and in the interstel- 
lar space are small. However, if the magnetic field in 
the PWN/SNR is strongly amplified from the interstellar 
value, the energy loss rate due to the synchrotron emis- 
sion may be faster than that due to the adiabatic expan- 
sion of the SNR, and even than that due to the decline 
rate of the escape energy. In such case, the confined elec- 
trons/positrons cannot escape the SNR later and only 
the first component iV ei0SC ,i would be emitted and ob- 
served at the Earth. Anyway, as either the strength of 
the magnetic field in the PWN/SNR nor its time evolu- 
tion is generally uncertain, it is difficult to estimate the 
radiative cooling rate with a moderate strength of the 
magnetic field in a reliable way. In the following, we in- 
vestigate only the case that the radiative energy loss is 
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not significant. 

4. OBSERVED ELECTRON SPECTRUM IN THE ESCAPE 
SCENARIO 

Now that we have the time evolution of the escap- 
ing CR particle flux from the source, we can obtain the 
observed electron spectrum by solving the propagation 
of CR electrons/positrons with the diffusion equation 
shown in Eq.Q. Once we know the Green's function 
of this equation with respect to the time and position, 
G(t, r, e e ; r), we can obtain the observed electron spec- 
trum as 



f(t,r,e e )= / G(t,r,e e ;T)dT, 



(20) 



where ti is the time when the particle injection has 
started, which is assumed to be equal to isedov in the 
following discussions. 

The mathematical description of G{t,r,e e ;r) was de- 
rived by Atoyan et al. (1995), 



G(t,r,s e ;t ) 



TT 3 / 2 P(e e )d d itf{e efi )'- 



■ exp 



ddiff (£e : o) 



m) 



where e e ,o is the energy of electrons/positrons at the time 
t which arc cooled down to e e at the time i, and d^is is 
the diffusion length given by 



rfdiff — 2 



! -° D(x)dx 



-1/2 



(22) 



as shown in Eqs.(10) and (11) in Atoyan et al. (1995). 

In deriving the energy loss rate P(e e ), we use the for- 
mulation shown by Moderski et al. (2005) 



3m 2 c 3 



8tt 



de 7 Mtot(e 7 )/KN 



4:£ e E>y 
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where ot is the Thomson cross section, u t ot(£ 7 )cfe 7 is 
the energy density of interstellar photons with the en- 
ergy between e 7 and e 7 + de 1 (including CMB, starlight 
and dust emission; Porter et al. 2008), and B is the in- 
terstellar magnetic field which we here set as 1/iG. Here 
the function /kn(i) is the correction factor to include 
the Klein- Nishina effect. According to Moderski et al. 
(2005), this function can be expressed as 



/kn(6) 
where b = 4£ e e 7 /(r7j e c 2 ) 2 
9(b) = 



9g(b) 



(24) 



2 



6 



6) 



^b 3 + 6b 2 + % + 4 



{1 + bf 



-2 + 2Li 2 (-6), 
and the function Li2(z) is the dilogarithm 

r0 ln(l - t)dt 



Li 2 (z) 



t 



(25) 



(26) 



According to the recent experiments, especially 
H.E.S.S. (Aharonian et al. 2008), the background CR 



electron/positron flux seems to have a high energy drop- 
ping around a few TeV. This dropping is quite naturally 
explained in the context of the astrophysical origin of CR 
electrons/positrons because the number of the sources 
contributing to the TeV energy band is quite small ac- 
cording to the birth rate of SNe/pulsars in the vicinity 
of the Earth (Kobayashi et al. 2004; Kawanaka et al. 
2010). In fact, since the pulsars which contribute to the 
electron flux at the energy e e should be younger than 
the cooling time of electrons/positrons t coo \ ~ l/(be e ) 
and should be located cl oser to the Earth than the dif- 
fusion length (idiff ~ 2-v/ K(e e )t coo \, the number of the 
pulsars contributing to > TeV band should be as small 



/ p \ -5/3 / 
N PSR (e e ) ~ 6 (^) ( 



R 



0.7 x lO-V^kpc" 



1(27) 



where R is the local pulsar birth rate per unit surface area 
of our Galaxy. If we can separate the contribution of a 
single young source from the observed electron spectrum, 
we can get the information of the CR injection into the 
ISM from that source. For this reason, we especially pay 
attention to the TeV spectral features of CR electrons 
from a pulsar in the followings. 

In Figs. 3 and 4 we show the time evolutions of CR 
electron spectrum from a nearby pulsar according to the 
models of the escape energy e esc (t) adopted in the previ- 
ous section (see also Fig. 2). It is clear that there exists a 
low energy cutoff in each spectrum corresponding to the 
value of e csc (t) at that time. The spectral shapes gener- 
ally depend on other parameters such as the high energy 
break of the intrinsic electron spectrum, the spectral in- 
dex, the duration of electron/positron injection from a 
pulsar, and the total energy of CR electrons/positrons. 
However, the sharp cutoff feature in the low energy side 
of the spectrum is almost independent of these proper- 
ties. In Fig. 3, we can see that the low energy cutoff of 
each spectrum is slightly broadened compared with that 
in Fig. 4. This is because the model adopted in this 
figure assumes that s esc (t) decreases more rapidly than 
in the case of Fig. 4 and so the CR electrons/positrons 
in the broader energy range can reach the observer while 
in the case of Fig. 4 where e csc (t) decreases slowly the 
low energy cutoff becomes very narrow. In either case 
the dropoff in the low energy side of the spectrum is so 
steep that one should assume the intrinsic spectral index 
as hard as a < — 1 if we neglect the energy-dependent 
CR escape effects. As we mentioned in the last section, 
if the magnetic field 

The CR electron spectrum with the age and distance 
similar to the Vela pulsar (i agc ~ 10 4 year, r ~ 290pc), 
which is thought to be surrounded by the supernova rem- 
nant (Aschenbach et al. 1995), is shown in Fig. 5. Here 
we show the spectrum with the escape model of Ptuskin 
& Zirakashvili (2005) as well as the spectrum without 
the CR confinement in the SNR. In addition, the elec- 
tron flux which have not been confined in the SNR (i.e. 
N e . 0S c,i) and that which have once confined and escaped 
later (N e .csc.2) are shown. We can see that the latter 
component dominates the flux around the low energy 
cutoff. In this case we can detect the electron flux from 
Vela pulsar with such a sharp spectral cutoff by near 
future missions such as CALET (we assume the geomet- 
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Fig. 3. — The electron spectra from a single pulsar surrounded 
by a supernova remnant with different ages; 3 X 10 3 (dotted line), 
5 X 10 3 (dashed line), 1 X 10 4 (solid line) and 2 X 10 4 years (dot- 
dashed line). As for the model of the escape threshold energy 
£esc(i) we adopt the results by Gabici et al. (2009; see Fig. 2). 
We adopt the background model of exponentially cutoff power-law 
with an index of -3.0 and a cutoff at 1.5TeV, which is similar to 
that shown in Aharonian et al. (2008) and reproduces the data 
in ~10GeV-lTeV well, and we show the spectrum including this 
background by thick lines. We assume that a source at r = 290pc 
from the Earth produces pairs with total energy E + = E e - = 
0.5 X 10 48 erg , duration to = 10 4 year, spectral index a = 2.0 and 
the high energy break £,.,cut = lOTeV. 
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Fig. 4. — The same plots with Fig. 3., but using the escape model 
by Ptuskin & Zirakashvili (2005; see Fig. 2). 
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Fig. 5. — The electron spectrum from a single pulsar surrounded 
by a supernova remnant (thick solid line). As for the model for 
the escape threshold energy e CS c(t) we adopt that of Ptuskin & 
Zirakashvili (2005; see Fig. 2), and the background model (thin 
dotted line) is the same as used in Fig. 3 and Fig. 4. The spec- 
trum without assuming the energy-dependent escape (thin dashed 
line), the flux from N e ,esc,l and iV e ,csc,2 (long-dashed line and 
dot-dashed line, respectively), and the errorbars expected from the 
5-years observation by CALET (5T2T = 220m 2 sr days; red down- 
ward triangles) are also shown. 



rical factor Sfl times the observation time T ~ 5years as 
~ 220m 2 sr days) because the assumed electron/positron 
flux is sufficiently large and the low energy cutoff comes 
beyond the high energy dropping of the background flux 
inferred by H.E.S.S. (Aharonian et al. 2009). If such 
a sharp cutoff feature is confirmed, that would be the 
strong evidence of the energy-dependent escape of CR 
particles from the acceleration site. As we mentioned 
in the last section the second component iV ej esc,2 would 
not be observed if the magnetic field in the PWN is 
strong and the electrons/positrons once confined in the 
PWN/SNR lose their energy too rapidly to escape the 
SNR later. However, the first component iV ei esc,i also 
has a sharp cutoff in the low energy side of its spectrum 
and so even in this case it is still possible to prove the 
CR escape from an SNR from the electron spectrum. If 
the threshold energy £esc(i) at that time is smaller than 
a few TeV then the cutoff feature would be hidden by the 
background flux. In fact, assuming the phenomenologi- 
cal model of e csc (t) adopted by Gabici et al. (2009) and 
Ohira et al. (2010a), the low energy end of the spectrum 
from the Vela pulsar would be buried in the background 
flux because the escape energy at the age of Vela SNR 
(~ 10 4 year) becomes smaller than ~ a few TeV, and so 
the cutoff feature due to the energy-dependent escape of 
CR electrons/positrons may not be resolved. 

As we mention in Sec. 1, the CR escape from a source 
(i.e. SNR) seems to be an important process for the ob- 
served CR spectrum below the knee (~ 10 15 5 eV) and 
the broken power-law gamma-ray spectra observed from 
SNRs interacting with molecular clouds. As for the spec- 
trum of CR protons and nuclei we can see only the su- 
perposed flux from multiple CR sources with different 
ages, and so it is impossible to see the direct evidence 
of energy-dependent CR escape from a single CR source. 
On the other hand, in the ~ 1 - lOTeV band of the CR 
electron spectrum we may be able to get the spectrum 
from a single source, and then we expect the sharp low 
energy cutoff showing the boundary between runaway 
CR electrons/positrons and confined CRs whose energy 
is not high enough to escape from the SNR shock. There- 
fore, this cutoff can be the first direct evidence that the 
energy-dependent CR escape really occurs at the source 
of those CR electrons/positrons. 

5. SUMMARY AND DISCUSSIONS 

In this paper we show the possibility of getting the 
evidence of the energy-dependent CR escape from the 
SNR shock by observing the TeV CR electron spectrum. 
The escape of CR particles is the fundamental process in 
emitting the CR particles accelerated at the SNR shock, 
and the history of the escape energy e eS c(t), which is the 
threshold energy of accelerated particles not to be con- 
fined but to escape the shock into the ISM, is essential in 
determining the observed CR spectrum. Although there 
are some indirect suggestions to this scenario from recent 
gamma-ray observations of SNRs interacting with molec- 
ular clouds, we have not verified the energy-dependent 
CR escape model of SNRs directly from the observations. 

In order to see the history of CR injections from SNRs 
we focus on the lepton component of CRs. In the en- 
ergy band larger than ~ TeV, it is expected that the 
electron/positron flux from a few nearby young astro- 
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physical sources (e.g. Vela pulsar) can be observed, then 
we will be able to get the information about the intrin- 
sic CR spectrum from a single source at a certain time, 
which has not been well understood because the observed 
nuclear CR spectrum consists of the contributions from 
multiple sources and it is impossible to resolve the spec- 
trum from a single source. Especially if a nearby young 
pulsar embedded in an SNR emits sufficiently large flux 
of CR electrons/positrons, they are injected into the SNR 
shock and only the electrons/positrons with energy larger 
than e esc (t) can go through the SNR and propagate into 
the ISM. As a result, the observed electron spectrum 
from that pulsar may have the low energy cutoff corre- 
sponding to the escape energy e esc (t) at the observation 
time. In order to detect such spectral feature the low 
energy cutoff (i.e. e esc (t) at the observation time) should 
come above > 1 — lOTeV because otherwise the cutoff 
would be buried in the background flux and would be 
hard to resolve. 

The possibility that there is such a low energy cut- 
off in TeV electron spectrum has firstly pointed out by 
this study, and taking into account the variation of the 
CR escape history from the SNR and/or the age of the 
pulsar, we can predict a variety of spectral shapes that 
have not been considered in the context of astrophysi- 
cal sources. For example, if the high energy cutoff of 
the spectrum, which is determined by the intrinsic cutoff 
energy at the pulsar and the radiative energy loss of elec- 
trons/positrons during their propagation, is nearly equal 
to the low energy cutoff determined by the escape from 
the SNR, the resulting spectrum would have a narrow- 
line like shape and it would be the similar feature to 
that expected from the annihilation of dark matter parti- 
cles. Whether the observed CR electron spectrum is orig- 
inated from a single nearby source would be determined 
by measuring an anisotropy of CR electrons/positrons 
(Mao & Shen 1972; Buesching et al. 2008; Ioka 2010). 

Note that the high energy break energy of the electron 
spectrum from a pulsar £ e , C ut generally depends on time 



Abdo, A. A. et al. 2009a, Phys. Rev. Lett. 102, 181101 
Abdo, A. A. et al. 2009b, ApJ, 706, LI 
Abdo, A. A. et al. 2010a, ApJ, 712, 459 
Abdo, A. A. et al. 2010b, ApJ, 718, 348 
Abdo, A. A. et al. 2010c, Science, 327, 1103 
Ackermann, M. et al. 2010, arXiv:1008.3999 
Adriani, O. et al. 2009, Nature, 458, 607 
Aharonian, F. A. et al. 2004, Nature, 432, 75 
Aharonian, F. A. et al. 2005, A&A, 437, L7 
Aharonian, F. A. et al. 2006, A&A, 448, L43 
Aharonian, F. A. et al. 2008, Phys. Rev. Lett. 101, 261104 
Aharonian, F. A. et al. 2009, A&A, 508, 561 
Aharonian, F. A. & Atoyan, A. 1996, A&A 309, 917 
Ahlers, M., Mertsch, P. & Sarkar, S. 2009, Phys. Rev. D 80, 
123017 

Aschenbach, B., Egger, R., & Trumper, J. 1995, Nature, 373, 587 
Atoyan, A. M., Aharonian, F. A. & Volk, H. J. 1995, Phys. Rev. 
D 52, 3265 

Beischer, B., von Doetinchem, P., Gast, H., Kirn, T. & Schael, S. 

2009, New J. of Phys. 11, 5021 
Biermann, P. L., Becker, J. K., Meli, A., Rhode, W., Seo, E. S. & 

Stanev, T. 2009, Phys. Rev. Lett. 103, 061101 
Blandford, R. & Eichler, D. 1987, Phys. Rep., 154, 1 
Blasi, P. 2009, Phys. Rev. Lett. 103, 051104 
Blasi, P. & Serpico, P. D. 2009, Phys. Rev. Lett. 103, 081103 
Blasi, P. & Amato, E. 2010, arXiv: 1007.4745 



and should be determined from the analysis of the par- 
ticle acceleration processes in the pulsar wind nebula as 
well as the radiative cooling of electrons/positrons dur- 
ing the acceleration. However, these processes have not 
been fully understood yet and so it is difficult to give the 
maximum electron/positron energy from the first princi- 
ple. For the present purpose, we are interested in the 
existence of the low energy cutoff due to the energy- 
dependent CR escape and unless the intrinsic high en- 
ergy break e e ,cut is much smaller than the escape energy 
e esc at the observation time, there would be a sufficiently 
large electron/positron flux above e csc and so the cutoff 
feature would be detected clearly enough to probe the 
CR escape scenario. 

We should also mention that the position of the low 
energy cutoff depends on the model of the time evolu- 
tion of the CR escape energy and the highly unknown 
magnetic and radiation fields, which produce the energy 
losses of electrons/positrons. 

The energy- dependent escape of CR elec- 
trons/positrons may be confirmed by the observations 
of the radio to 7-ray emissions from an SNR with a 
pulsar. The electrons/positrons escaping the SNR would 
radiate via synchrotron emissions and inverse Compton 
scatterings just outside the SNR shock. If we can 
observe the radio to 7-ray intensity around the nearby 
young SNR with a pulsar and can fit that spectrum by 
the leptonic model with the electrons/positrons having 
a low energy cutoff in their energy distributions, that 
would be the support of the energy-dependent escape of 
CRs from the SNR (in preparation). 
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